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ABSTRACT 

The three-point correlation function (3PCF) of the spin-2 fields, cosmic shear and microwave background 
polarization, is a statistical measure of non-Gaussian signals. At each vertex of a triangle, the shear field has two 
independent components. The resulting eight possible 3PCFs were recently investigated by Schneider & Lombardi 
(2002) and Zaldarriaga & Scoccimarro (2002). Using rotation and parity transformations they posed the question: 
how many components of the shear 3PCF are non-zero and useful? We address this question using an analytical 
model and measurements from ray-tracing simulations. We show that all the eight 3PCFs are generally non-zero 
and have comparable amplitude. These eight 3PCFs can be used to improve the signal-to-noise from survey data 
and their configuration dependence can be used to separate the contribution from E- and B-modes. This separation 
provides a new and precise probe of systematic errors. We estimate the signal-to-noise for measuring the shear 
3PCF from weak lensing surveys using simulated maps that include the noise due to intrinsic ellipticities. A deep 
lensing survey with area of order 10 square degrees would allow for the detection of the shear 3PCFs; a survey 
with area exceeding 100 square degrees is needed for accurate measurements. 

Subject headings: cosmology: theory - gravitational lensing - large-scale structure of universe - cosmic 
microwave background 



1. INTRODUCTION 

Can the three-point correlation function (3PCF) of a spin- 
2 field be a useful probe of non-Gaussian signals in cosmol- 
ogy? Most previous work has focused on the 3PCF or its 
Fourier transform, the bispectrum, of scalar quantities such as 
the galaxy number density and temperature anisotropies in the 
cosmic microwave background (CMB). An exciting develop- 
ment is the recent detection of the 3PCF of the cosmic shear 
field by Bernardeau, Mellier & Van Waerbeke (2002a; also see 
Bernardeau, Van Waerbeke & Mellier 2002b), which provides 
new constraints on structure formation models beyond those 
provided by two-point statistics. For the CMB, the bispectrum 
of the temperature anisotropies is a promising way of probing 
primordial non-Gaussian perturbations (e.g., Verde et al. 2000; 
Komatsu et al. 2002). The 3PCF of the CMB polarization 
might open a new, complementary window for this in that it 
can separately measure the non-Gaussian signals arising from 
primordial scalar, vector and tensor perturbations. 

A spin-2 field has two components on the sky, hence its 3PCF 
in general has 2 3 = 8 components. One may ask: how many of 
these components are non-zero? Do all components carry use- 
ful information about the B- and Zi-modes of the spin-2 field 1 ? 
Recently, Schneider & Lombardi (2002; hereafter SL02) and 
Zaldarriaga & Scoccimarro (2002; hereafter ZS02) investigated 
these questions. The purpose of this Letter is to clarify these is- 
sues using geometrical arguments and measurements from ray- 
tracing simulations of the cosmic shear. We will pay particular 
attention to the problem of how the eight 3PCFs are related to 
the E- and Z?-modes. 

2. THE 3PCF OF SPIN-2 FIELDS 

The two components of a spin-2 field depend on the choice 
of the coordinate system. Suppose we have the shear compo- 

1 A two-dimensional spin-2 field can be separated into an £-mode derivable 
from a scalar potential and a pseudoscalar B-mode (Kamionkowski, Kosowski 
& Stebbins 1997; Zaldarriaga & Seljak 1997; Hu & White 1997 for the CMB 
polarization and Stebbins 1996; Kamionkowski et al. 1998; Crittenden et al. 
2001; Schneider et al. 2002 for the cosmic shear). 



nents, 71 and 72 on the sky, for given Cartesian coordinates . 
A rotation of the coordinate system by ip (in the anticlock- 
wise direction in our convention) transforms the shear fields as 

For the two-point correlation function (2PCF), the prob- 
lem of the coordinate dependence of the shear field has been 
well studied in the literature (e.g., Bartelmann & Schneider 
2001). For a given pair of points, Xi and X2, separated 
by a fixed angle 6, we can define two components of the 
2PCF which are invariant under coordinate rotations: £+(9) = 
( 7+ (X 1 ) 7+ (X 2 )) and £ x (0) = (7x(X I ) 7x (X 2 )). Here 7+ (X,-) 
and 7><(X,) are the shear components defined by projecting 
the shear along the direction \p connecting the two points: 
7+ (X,) + i lx (X;) = - [71 (X;) + i 72 (Xi)] e- 1 ^. The other possible 
component, (7+7x), vanishes because of parity invariance. 

SL02 and ZS02 investigated possible eight components of 
the shear 3PCF based on the +/x decomposition of the shear 
field. We will closely follow these authors in setting up the ge- 
ometry of the 3PCF. In contrast to the 2PCF, there is no unique 
choice of a reference direction to define the +/x decomposi- 
tions for the shear fields at the vertices of a triangle. The 'cen- 
ter of mass', o, of a triangle is one possible fiducial choice: 
o = (1 /3)J^Li X,-. Figure 1 shows a sketch of the triangle we 
use to define the shear 3PCFs. The solid and dashed lines at 
each vertex show the positive directions of the + and x compo- 
nents, respectively. We also define the interior angle ip between 
xi and X2. The projection operator to compute the + or x com- 
ponents at each vertex is written as: 

P + (Xd = -(0fi-0f 2 ^a)/0f, 

p x (x ( ) = -(-2%^, e\ - e\)/el (i) 

where <?,■ = X, - o. Using these projections, we obtain at each 
vertex 7+ = P +l jx +.F+272 and 7* = P x 171 +P X 2j2- These trans- 
form under parity as 7+ — > 7+ and 7 X — *■ — 7 X . In terms of 7+ 
and 7 X we define the eight components of the shear 3PCF for a 

2 We use the notations of weak lensing for the spin-2 field, 7,; our discussion 
can be applied to the CMB polarization if one replaces 71 and 72 with the 
Stokes parameters Q and U, respectively. 
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FIG. 1 . — Upper. Definitions of the triangle variables used for the 3PCF. The 
interior angle between x 1 and x 2 is defined as ip and its positive direction is 
shown by the arrow. Note that the solid and dashed curves at each vertex show 
the positive directions for 7+ and -y x . Lower: The mirror transformation with 
respect to the side vector Xi for f++ x . The vertices are transformed as 1 — ► 1' 
and so on. Note that the + components remain unchanged in sign, but the sign 
of the X component is flipped. 



given triangle as a function of x\, x 2 and ip: 

W*i,*2,V) = (7p(Xi) 7 ,(X 2 )7r(X 3 )), (2) 

where /1, is, t = + or x and (• • •) denotes the ensemble average. 
The 3PCF defined above is invariant under triangle rotations 
with respect to the center o as pointed out in SL02. Hence it 
is fully characterized by the three variables x\,x 2 ,ip, just like 
the 3PCF of a scalar field (note that if we had used the Carte- 
sian components of the shear, we would need four variables to 
specify it as it is not invariant under rotation). 

For a pure E-mode, based on the properties of 7+ and j x un- 
der parity transformations described above, we divide the eight 
shear 3PCFs into two groups: 

Parity-even functions: C+++, C+x x , Cx+x , Cx x+ 

Parity-odd functions: Cxxx,Cx++,C+x+,C++x- (3) 

Weak lensing only produces the E-mode, while source galaxy 
clustering, intrinsic alignments and observational systemat- 
ics induce both E and B-modes in general (Crittenden et al. 
2001; Schneider et al. 2002) 3 . For the CMB polarization, 
although primordial scalar perturbations generate only the E- 
mode, vector and tensor perturbations can induce both modes 
(e.g., Kamionkowski et al. 1997). 

A question posed by SL02 and ZS02 is whether all eight 
3PCFs could have contributions from lensing. We examine this 
question with geometrical considerations and show results from 
simulations in the next section. The parity transformation for a 
triangle can be taken to be a mirror transformation as sketched 
in Figure 1, which shows the parity transformation with re- 
spect to the side vector Xi as one example. It corresponds to 
the change tp — > 2ir-tp in our parameters. Figure 1 illustrates 
the transformation for C++x- From statistical homogeneity and 
symmetry, the amplitude of the 3PCF depends only on the dis- 
tances between the center and each vertex. Hence, the absolute 
amplitudes of (++ x for the two triangles shown should be the 
same. But the sign of j x at the vertex 3' changes under this 

'Multiple lensing deflections generally induce the B-mode, but Jain, Seljak & 
White (2000) showed that the induced B-mode is much smaller than the lensing 
£-mode using ray-tracing simulations. 



parity transformation. For the eight 3PCFs in equation (3), we 
can say in general: 

Parity-even: (xi,x 2 ,ip) = C^r (xi,x 2 ,2w-ip), 

Parity-odd: Cawt(*i,*2, if>) = ~C^r (x\,x 2 , 2n~i>). (4) 

Next, we consider an isosceles triangle with x\ = x 2 . In this 
case, the j x at vertex 3 and 3' in Figure 1 are statistically iden- 
tical (viewed from the center of the triangle, they should have 
equal contributions when averaged over a matter distribution). 
We thus have the additional symmetries for the two parity-odd 
3PCFs: 

C++ X (x u x u ip) = C++ x 27T-V0, 
Cxxx(*l,*l,^) = Cx xx (Xi,*i,27T-V>). (5) 
The properties described in equations (4) and (5) lead to 

Isosceles Triangles: (++ x =0, ( xxx =0. (6) 

Note that this argument does not lead to the vanishing of the 
other two parity-odd functions, in which the x component is 
at a vertex bounded by unequal sides. For equilateral triangles 
however all four parity-odd functions vanish: 

Equilateral Triangles : Cparity-odd = 0. (7) 

We argue that for a generic triangle with all sides unequal, the 
differences in side lengths break the symmetry that is expressed 
in equation (5). In other words, the fact that gravitation has a 
scale dependence allows parity-odd 3PCFs to be non-zero for 
a triangle with unequal sides. These conclusions are consistent 
with those of SL02. 

3. PREDICTIONS FROM RAY TRACING SIMULATIONS 

To demonstrate the properties of the shear 3PCFs more pre- 
cisely, we employ ray-tracing simulations of the cosmic shear 
performed by Jain et al. (2000). We use the SCDM model 
(fi m o =1, h = 0.5 and as = 0.6), with source redshift z s = 1 
and area fi s k y = 7.68 degree 2 (see Jain et al. 2000 for more 
details). We followed Barriga & Gaztanaga (2002) for the al- 
gorithm to calculate the 3PCF. Lists of neighbors are used to 
find the three vertices from the cell-based data as well as to cal- 
culate the projection operators for each vertex. The error bars 
shown in the following figures are computed from 9 different 
realizations (see Takada & Jain 2002c for more details). 

Figure 2 shows the results for the shear 3PCFs for isosce- 
les triangles against the interior angle ip between the two equal 
sides x\ = x 2 = 0'.97 (see Figure 1). The shear 3PCF is com- 
puted from the simulated data by averaging the estimator over 
all triplets with given triangle configuration. For the bin size we 
used Ax = 0.'08 for the side lengths and Aip = 7r/20 for the inte- 
rior angle. The upper panel shows the four parity-even 3PCFs, 
while the lower panel shows the parity-odd functions. As ex- 
pected, £+++ appears to have the greatest lensing contributions, 
and, in particular, peaks for equilateral triangles with ip = ir/3. 
The other three modes have smaller amplitude than £+++, and 
two of them are the same because of symmetric triangle and 
shear configurations. Interestingly, the features of these curves 
match the theoretical curves shown in Figure 3 of ZS02 (</> in 
their figures corresponds to ir-tp in our figure). Their results 
were presented with arbitrary normalization of the y-axis, so 
we can only compare the configuration dependence. Thus these 
small-scale complex features are basically captured by the tan- 
gential shear patterns around a single halo. 

From the lower panel in Figure 2, it is clear that two of the 
parity-odd 3PCFs are non-zero, but the others vanish, as stated 
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FIG. 2. — Ray-tracing simulation results of the shear 3PCF for isosceles 
triangles. The two side lengths are x\ = x% = 0'97. The upper panel shows the 
four parity-even functions, while the lower panel shows the four parity-odd 
functions. For illustrative purpose, the results for £ x +x and C++x are slightly 
shifted horizontally. The solid curves show analytical calculations of the shear 
3PCFs based on the halo model developed in Takada & Jain (2002b). 



in equation (6). This result clarifies that the parity-odd 3PCFs in 
general do carry lensing information and vanish only for special 
triangle configurations. For isosceles triangles, the two vanish- 
ing parity-odd functions can be used to discriminate systemat- 
ics from the E'-mode. For ip = ir/3, the triangle is equilateral, 
and all the parity-odd 3PCFs vanish as stated in equation (7). 

We have verified that the simulation results are in agreement 
with the 1-halo term predictions of the halo model. The solid 
curves in Figure 2 show the halo model predictions that fol- 
low the real space halo formalism developed in Takada & Jain 
(2002b). The calculations are based on adapting equation (52) 
from this paper by replacing K m with the relevant shear compo- 
nent for a halo of given mass. A detailed comparison of analyti- 
cal and simulation results for different models will be presented 
in Takada & Jain (2002c). 

The features shown above for isosceles triangles are explored 
for a triangle with three unequal sides in Figure 3, in which 
x\ =X2/2 = 0.'97. This figure is plotted over the range ip = [0,27r] 
to explicitly show the parity-even (upper panel) and odd (lower 
panel) properties of the shear 3PCFs. Note that the parity trans- 
formation is equivalent to the change ip — ► 27T — tp in our pa- 
rameters. As noted in ZS02, the amplitude of the 3PCF is sup- 
pressed for elongated triangles because of cancellations in the 
signals by the vector-like property of the shear. The figure re- 
veals that all the parity-odd 3PCFs are non-zero and in fact all 
the eight 3PCFs have roughly comparable amplitude. Note that 
Cx x x and (+x+ vanish at ip/w = 0.42, as the triangle becomes 
an isosceles triangle with xi = x$, and therefore the 3PCFs are 
consistent with the result in Figure 2. 

By rotating the simulated shear field at each position by 45 
degrees, we can investigate the 3PCF of a pure B-mode map 
(Kaiser 1992). Since this procedure transforms 7+ — > j x and 
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FIG. 3. — The results for a general triangle with x\ = X2/2 = 0f97, as in 
Figure 2. 
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FIG. 4. — The upper panel shows the shear 3PCF for equilateral triangle as a 
function of the side length. Ray tracing measurements and halo model calcula- 
tions are shown for the two independent, non-zero 3PCFs, (+++ and Cx++- The 
lower panel shows £+++ and £x x x from simulated maps that include noise 
due to intrinsic ellipticities with rms u f = 0.3. The bin size is Ar = 0'.325. 



7x we find -> Cx x x and so on - Hence, Cx x x 

carries most of the fi-mode signal for triangles that are close 
to equilateral. Further the symmetry properties of the 3PCFs 
under ip — > 2ir— ip (shown in Figure 3) are reversed for the B- 
mode; for example, Cx x x (VO = Cx x x - ip)- Therefore, we 
should keep in mind that for a general spin-2 field which in- 
cludes both E and B modes, the property ((ip) = ±C(27r - ip) 
can no longer be expected. The measurements of C thus has 
to be done for the full range ip = [0,27r] (unlike the case for 
the 3PCF of a scalar quantity like the density field or the CMB 
temperature fluctuations field). 

The upper panel of Figure 4 shows the results for the shear 
3PCF for equilateral triangles against the side length. For the 
measurements from simulated data, we used triangles with each 
side length in the range [r— Ar, r+Ar], with bin size Ar = 0'.325 
kept fixed. Since all the parity-odd 3PCFs vanish for an equilat- 
eral triangle, the results for C+++ and C+x x (= Cx+x = C+x x ) are 
shown. We also show the halo model predictions for the 1-halo 
terms of C+++ and C+x x by the solid and dashed curves, respec- 
tively. It is again apparent that the halo model predictions are 
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in good agreement with the simulation results. 

To estimate the signal-to-noise (S/N) for measuring the shear 
3PCFs from an actual survey, we need to account for the noise 
arising from the large intrinsic ellipticities of source galax- 
ies. The lower panel of Figure 4 shows the 3PCFs measured 
from simulated shear maps including the noise contamination. 
We assume that the intrinsic ellipticities have random orienta- 
tions with a Gaussian distribution for the amplitudes with rms 
a e = 0.3. The galaxies are taken to be randomly distributed with 
number density n ga \ « 38 arcmin" 2 . One can see that ( +++ can 
be detected from a survey with area of order 10 square degrees, 
provided the errors are dominated by statistical errors. The 
other 3PCFs are more difficult to measure for equilateral trian- 
gles. Here we also show Cx x x which is expected to be zero for 
an E-mode signal, so the error bars on it show the accuracy with 
which the B-mode can be constrained. These results in part ver- 
ify the detection of the shear 3PCF by Bernardeau et al. (2002a) 
from the cosmic shear survey of Van Waerbeke et al. (2001), 
although they used a different estimator of the shear 3PCF. 
The S/N for one particular triangle configuration with size 

r scales roughly as S/N cx ^ 7 (cr e / v /7j^i')" 3 r2s l ( r 2 ve yr 2 (Alnr) 3//2 , 
where Sl survey is the survey area. By combining information 
from triangles with different configurations and from the dif- 
ferent 3PCFs, we can improve the S/N (however, to interpret 
the measured 3PCFs one must estimate how correlated neigh- 
boring bins are, as discussed below). Note further that we have 
used the SCDM model with C+++ = 2-4 x 10" 6 for r = 1-5'; 
the S/N for the concordance ACDM model is lower because 
the signal is smaller, e.g. £+++ = 4-8 x 10" 7 (Takada & Jain 
2002c). This rough analysis leads us to conclude that future 
surveys with area > 100 degree 2 should provide measurements 
of the 3PCFs with S/N well over 10 in each bin if systematic 
errors are eliminated. 

4. DISCUSSION 

We have investigated the 3PCF of the cosmic shear field. 
We measured the eight components of the shear 3PCF from 
ray-tracing simulations and checked our results using analyti- 
cal calculations based on the halo model. Figures 2 and 3 show 
that in general all eight components are non-zero, have com- 
parable amplitude, and have a complex configuration depen- 
dence. These results verify the analysis of SL02 who analyzed 
the 3PCFs based on their transformations under parity and ro- 
tations. 

The 3PCFs show non-Gaussian signatures induced by non- 
linear gravitational clustering. In constraining models from 
measurements of the shear 3PCFs, it will be necessary to de- 
velop an optimal way to combine the eight 3PCFs. To do this, 
we will need explicit relations showing how the shear 3PCFs 
are related to the lensing E-mode, in analogy with two-point 
statistics (Crittenden et al. 2001; Schneider et al. 1998, 2002). 
This problem is more easily formulated in Fourier space. If we 
use the Cartesian components of the shear to define the 3PCFs, 
then simple relations exist between the bispectra of the shear 
3PCFs, Qijk and the bispectrum of the convergence, B K : 

Cijk(h,h,h)=B K (luh,h)u i (h)u j (h)u k (h), (8) 

where i,j,k = 1,2 and m,(1) = (cos2<,oj, sin2cpj) for the E-mode. 
This equation shows that the shear 3PCFs can be related to a 
single 3PCF for the E-mode field. For example we can obtain 
the convergence bispectrum using the eight measured shear bis- 



pectra as: 

Bk(1i,12,1 3 ) = ^mK1i)m/12)^(1 3 )Cvi(1i,12,13). (9) 

i,j,k 

However in real space it is still unclear how best to use the eight 
3PCFs measured over a limited range of scales to get an opti- 
mal estimator of the E-mode 3PCF. Another question we have 
not considered is how strongly the shear 3PCFs are correlated 
for different triangle configurations, since lensing fields at ar- 
cminute scales are highly non-Gaussian (e.g. Takada & Jain 
2002a). 

Realistic data has noise which includes both E/B modes, as 
would contributions from intrinsic ellipticity correlations, non- 
linear lensing effects and systematic errors. We have shown 
that two or all four parity-odd 3PCFs vanish for isosceles or 
equilateral triangles, respectively (see Figures 2 and 4). For 
general triangles, Figure 3 shows that the 3PCFs of E-mode 
shear fields have specific symmetry properties under the change 
ip — > 2ir - ip, where ip is the interior angle of the triangle 
shown in Figure 1. These properties can be used to find B- 
mode contributions (whose symmetry properties are reversed) 
as a function of scale and configuration. Thus the origin of 
the B-mode contribution can be identified more precisely than 
is possible with two-point statistics. An example of an ex- 
plicit test for B modes is to measure combinations such as 
( ++x (xi,X2,tp) + C++x(xi,X2,2tt - ip). This should be zero for 
all the functions that are parity-odd for a pure E mode, as given 
in equation (3). 

The results presented above can also be applied to the Q and 
U Stokes parameters of the CMB polarization. The results we 
have shown imply that the eight 3PCFs constructed from com- 
binations of Q/U, e.g. (UQQ), have non-vanishing signals if 
the CMB polarization field is non-Gaussian. The 3PCFs have 
the advantage that it can discriminate the non-Gaussianity aris- 
ing from the E- and Z?-modes. However, since the polarization 
signal is small compared to the CMB temperature fluctuations, 
this will be a great challenge. 

We are grateful to G. Bernstein, L. Hui, M. Jarvis, P. Schnei- 
der, R. Scoccimarro, A. Stebbins, A. Szalay and M. Tegmark 
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NAG5-10923, NAG5-10924 and a Keck foundation grant. 
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